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Traffic  Modelling 


Activity  Types 


Commercial  shipping 
Commercial  fishing 
Commercial  recreational 
Aquaculture 
Ferries 


Cruise  ships 
Private  recreational 


Mill,  Jr.  W9BFtf 

A 

TvmF^EiDam  tmiditv  1  1mFF$&\  I  |- 

Class 

Available  Information 

Example 

A 

Continuous  time-stamped  position  data; 
vessel  info 

Canadian  West  Coast 
shipping  transits;  Automatic 
Identification  System  (AIS) 

B 

Specific  origin  &  destination,  and 
intermediate  waypoints;  travel  dates; 
vessel  info 

Canadian  East  Coast  shipping 
transits 

C 

Specific  origin  &  destination; 
travel  dates;  vessel  info 

Ferries;  Cruise  ships 

D 

Specific  origin;  general  destination; 
travel  dates;  vessel  info 

Commercial  fishing  by 

NAFO  zones 

E 

Specific  origin;  general  destination; 
frequency  of  trips;  vessel  info 

Lobster  fishing;  ecotours 

F 

General  origin;  general  destination; 
frequency  of  trips;  classes  of  vessels 

Recreational  boating 

Traffic  Modelling 


•  Maximize  use  of  information: 

b  source:  location,  date,  time 
a  intermediate  points:  way  points,  fishing  efforts 
b  final  destination:  location,  date,  time 

•  Extrapolate  using  more  general  information 

a  feasible  range:  distance,  time 

a  feasible  location;  eg.  fishing  grounds;  tourist 
sites 

a  typical  locations:  cluster  analysis 
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Dispersion  Algorithm 


Trips  to  feasible  region 


Grid  the  area 


O  Define  feasible  destinations 


generate  tracks 


Grid  results 


Abnormal  ship  trajectories 


Objectives 


(1)  To  mathematically  define  commonly  used 
shipping  routes  from  historical  shipping 
data. 

(2)  To  detect  unusual  vessel  trajectories 
based  on  norms  of  typical  shipping  routes 
using  Statistical  Process  Control  (SPC). 
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I  UIBHf  I  J  ^PJ  I  I  s4  .rgyi 

pattern  definitions 


©  Line  pattern  analysis  in  this  research  involves  the 
description  of  the  relationship  among  different 
trajectories  and  the  classification  of  similar 

trajectories  into  routes  and  paths. 

©  Trajectory:  a  line  generated  from  a  series  of 
waypoints  by  an  individual  moving  object. 

©  Route:  a  frequently  used  track  followed  by  moving 
objects  from  the  time  when  they  enter  the  study  area 
until  the  time  they  exit. 

©  Path:  a  feature  grouping  common  sections  of  routes. 
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Traffic  Monitoring  System 
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Geometric  Algorithm  Review 


Three  main  steps: 

•  Route-trajectory  matching; 

•  Route-trajectory  updating; 

•  Route- route  merging. 

®  Each  route  node  here  is  characterized  by  a 
vector  format  with  coordinates  x,  =  [x, ,  v, ]. 

©  Attributes  of  each  route:  central  axis,  left  and 
right  boundaries  (also  called  "envelopes"), 
normal  vector  and  weight. 

•  Resample  distance  -  equally  spaced  nodes  are 
obtained  on  each  trajectory  and  route. 
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Route-Trajectory  Matching 


•  Label  each  trajectory  with  a  unique  identifier  and 
resample  it.  The  first  trajectory  is  initialized  as  the  first 
existing  route  in  the  route  database.  . 

\ 

\ 

\ 


\ 

\ 


\ 

\ 

\ 

\ 
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Route-Trajectory  Matching 


•  A  second  trajectory  is  taken  from  the  database, 

resampled,  and  Hausdorff  distance  is  calculated  between 

the  trajectory  (T)  and  all  the  existing  rputes  (R). 

\ 

\ 

\ 

\ 
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Route-Trajectory  Matching 


9 


Compare  the  distances  d|(TR)  from  the  trajectory  to  all 
the  routes  in  the  route  database  and  get  the  minimum 

distance.  The  index  of  this  route  is  written  as  I . 

\ 

\ 
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Route-Trajectory  Matching 


®  I  f  d,  (TR)  is  greater  than  a  given  threshold 
( 1  trajectory  -  assimilate-  threshold' ) ,  the 
trajectory  will  be  created  as  a  new  route. 

•  Otherwise,  a  match  between  the  trajectory  and 
route  I  is  detected.  Go  to  2nd  step,  route 
updating. 
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Rout^T rajectory  Updating 


•  Central  axis  node: 


XIk:= 


w, 


Wj  +  1 


xIk  + 


1 


Wj  +  1 


X 


Ikt 


l> 


\ 
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End  Normal 
Line 
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Route-Jrajectovy  Updating 

•  Route  extension:  The  longer  part  of  the  trajectory  is  used 
to  extend  the  route  I .  I  f  the  extension  length  is  more  than 
half  of  the  length  of  the  entire  route,  the  trajectory  is 

created  as  a  new  route.  Extended  Nodes 


Rpute-Trajectory  Updating 

Boundaries:  The  farthest  intersections  on  both  sides 
are  the  new  left  and  right  boundary  nodes. 


•  Weight  of  route:  incremented  by  1. 


Extended  Nodes 


Route-Route  Merging 

•  Sort  all  the  routes  by  descending  order  of  weights. 

•  Step  1:  Start  from  the  route  with  the  highest  weight  (F^,). 
Calculate  the  Hausdorff  distance  between  Rh  and  the  route 
with  the  smallest  weight  in  the  database  (Rg). 

•  Step  2:  If  the  distance  is  less  than  a  predefined  threshold 
( ' route- merge- threshold' ) .  then  the  secondary  route  is 
qualified  to  be  merged  into  the  main  route  R^  Otherwise, 
they  are  not  merged. 

•  Step  3:  Go  to  the  next  lower  weighted  route  and  repeat  the 
above  two  steps  until  all  the  routes  in  the  database  are 
compared  with  R^ 

•  Step  4:  Go  to  the  next  higher  weighted  route  and  start  a 
new  loop  (the  above  3  steps),  until  all  the  routes  (except 
the  least  weighted  route)  are  used  as  Rh  to  be  compared. 


Wker  filexontrol 


'*  Route  Formation 


Path  Data 


Constants 


Run 


Output  Path  Data:  |D:\research\Route  Formation\RouteFornnation!2008-06-04_1 547]\PathsDutf2008-06-09_21 56). si  Browse 

Output  Junction  Data:  Browse 


T rajectory  Data 

1“  Skip  T rajectory  Data  Loading  (Route  Merge  and  Path  Formation  Only) 

Input  Trajectory  Data: 

(*  Text  C  Shape  Browse  | 

Trajectory  ID:  I - 3  F  Res-ample  Trajectories 

Output  Trajectory  Data:  |  D:\research\R  oute  Formation\RouteFormation|2003-0G-04_1 5471\Traput(2003-06-09_21  SGl.shp  Browse  | 

7  Calculate  Matches  for  T rajectories 

Output  Statistics  Data:  | D : \research\R oute  Formation\R outeFormationf 2008-06-04_1 547]\M atchS tatistics(2008-0G-0S_2'  _  Browse  | 

Test  Trajectories  (Text):  Browse  | 

Flagged  Test  Trajectories:  [da  —  Browse  | 

Route  Data 

17  Use  Input  Route  Data:  |  D:\research\R  oute  Formation^  utput[2008-06-06_1 1 22]_S orting\R outesO utf2006-0G-06_1 1 22l_  Browse  | 

Output  Route  Data:  [d!^  Browse 
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Constants 


User  parameter  control  (metres) 


ENVELOPE_OFFSET ■ 


M I N_T RAJ E CT 0 RY_LE NGTH_IN_RE SAM PLE  S  - 
P_Q_INTERSECTION_DIST_THRESHOLD  ■ 
MIN_C0MM0N_PR0P0RTI0N_F0R_MERGING  ■ 
MAX_UNC0MM0N_PR0P0RTI0N_F0R_MERGING  - 

C^^ECTO  R  Y_AS  S I M I LAT  E_T  H  RESH0LD> 
ROUTE_TRAJ_UPDATE_THRESHOLD  ■ 
ROUTE_TRAJ_BOUND_THRESHOLD  ■ 


foT" 

[2000 

r- 

[3000 

[O50” 

[Q50” 


[2500 

[2500 

[2500 


CJOUT  E_M  E  R  G  E_T  H  R  E  SH0LD> 
ROUTE_ROUTE_UPDATE_THRESHOLD 
ROUTE_ROUTE_BOUND_THRESHOLD 

MAX_ROUTE_EXTENSION_WEIGHT  ■ 

TRIM_N0DE_PR0P0RTI0N  ■ 
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Implementation  of  the  algorithm 


♦  Data  source:  VTOSS  database  (Vessel 
T raffic  Operations  and  Support  Systems) 


O  Study  area:  the  ships  that  approach  within 
50km  of  the  British  Columbia  (BC)  shoreline. 


♦  Study  period:  Fourteen  daily  datasets  from 
J  an.  1st  to  14th,  2003;  three  weekly  datasets 
-  J  an.  1st  to  7th,  8th  to  14th,  and  Aug.  6th  to 
12th,  2003. 


♦  The  following  figures  are  displayed  for  a 
one-week  traffic  data  (from  J  an.  1st  to  7th, 
2003) . 


Historical  trajectories 


(resample -interval  =  2500,  trajectory _assimilate -threshold 
=  3000,  route jnerge -threshold  =  3000) 


Kernel  Density  of  one  week  of  traffic 


lanaimo 


Victor! 


SPC  Methodology  Review 

•  Cumulative  sum  (CUSUM)  control  chart  is 
applied  in  maritime  traffic  to  monitor  the  process 
(i.e.  movement)  of  any  new  individual  trajectory. 

•  Let  x,  be  the  ith  observation  on  the  process.  |u0  is 
the  target  value  for  the  quality  characteristic  x. 

•  The  tabular  CUSUM  works  by  accumulating 
deviations  from  |u0  that  are  above  target  with 
one  statistic  C+  and  accumulating  deviations 
from  |u0  that  are  below  target  with  C. 

C-  =  max[0,  x,  -  ( ju0  +  K)  +  C+,  ] 

C:  =  max[0,  (ji0  -K)-  x,  +  C:,  ] 
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^j^f^Iethodqlogy  Revie i| 


►  Xj  =  the  minimum  distance  from  node  /on 
the  trajectory  to  the  route  envelope. 

&The  overall  mean  and  standard  deviation 
(unbiased)  are  calculated  for  the  route 
across  all  incorporated  trajectory  nodes  to 
route  envelope  distances. 

>  Only  C+  statistic  is  used,  because  distance 
is  a  positive  variable. 

*Two  levels  of  alarm  are  set. 
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individual  trajectories 


Traj:  5-380 


Traj:  2-516 


Discrim  in  ating 
recreational  boat  types 


my  |  |  '■■or  I  I  I  i 

Characterizing  of  Recreational  Boats 


based  on  GPS  Trajectory  Points 
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Objectives 


This  Study  examines  Global  Positioning  System 
(GPS)  trajectory  data  to  characterize  spatial 
patterns,  providing  insight  into  recreational  boat 
movements  and  amount  of  exposure  to  risk. 


Spatial  Patternrf 


>  Patterns  of  Different  Recreational  Boat  Types 


(a)  c-0707-08-03 


(b)  c-0721-11-04 


(c)  s-0722-14-10 


Spatial  Pattern  II 


>  Patterns  of  Different  Recreational  Boat 
Types  at  Different  Locations 


(a)  s-0722-14-10 


(b)  s-0824-16-03 


Scope  of  Study 


Sample  size 

NB 

NS 

Canoe 

4 

17 

Kayak 

8 

21 

Motorboat 

15 

10 

Sailboat 

5 

47 

Methodology 


>  Data  Pre-cleaning 

>  Points  on  land 

>  Time  gap 

>  Resting/Stopping  events 

>  Attributes  Extracting/Calculation 

> Speed 

>  Total  distance 

>  Bounding  Box 
>Dedensified  Trajectory 
>Distance  from  shore 

>  Discrimination  &  Classification 


Attributes  to  Variates 


> Speed 

>Max  Speed 

>Max  1/20  Speed 

>  Total  distance 

>Mean  Speed 

>  Bounding  Box 

- > 

>  Total  Distance 

K 

>  Complexity 

>Dedensified  Trajectory 

>  Aspect  Ratio 

>Distance  from  shore 

>Mean  Turning  Angle 

>  Distance  to  Shore 

Pattern  Analysis 


>  Discrimination 


>  Classification 


Discriminant  Analysis 

Do  different  boat  types  have  significant  different  attributes? 


Attribute 

NB 

NS 

TD 

CO  cPO 

Mean.S 

CO 

CCDCDCD 

Max.S 

CO 

CDDCD 

Max  1/20.  S 

CO 

CDDCD 

TAQ 

<?D>  (IT) 

CD  CD 

TAH 

cgiO  CO 

D>  CD 

C 

CO  CO 

DTS 

(C,K,M)  CO 

CD  CD 

AR 

CD  CD 

CD  CD 

Discriminant  Analysis 


Do  geography  characteristics  influence  boating  trajectories? 


Canoe 

Kayak 

Motorboat 

Sailboat 

TAH 

X 

X 

X 

TAQ 

X 

X 

X 

X 

c 

X 

X 

X 

X 

DTS 

X 

X 

X 

Max.S 

X 

X 

Max  i/2o* S 

X 

X 

Mean.S 

X 

X 

X 

TD 

X 

X 

X 

AR 

X 

X 

X 

P=0.1 

X:  Not  reject  H0:  There  is  no  significant  difference  between  the  two  areas. 


>  Problem:  Whether  can  we  determine  the  type  of  recreational 
boat  solely  through  the  trajectory  data  without  reference  to 
the  weather  and  operator  factors,  even  geography? 

>  Objectives: 

>  Determining  whether  statistically  significant  differences 
exist  between  attributes  of  different  boating  trajectories  of 
different  boat  types 

>  Determining  which  of  the  independent  variables  account 
the  most  for  the  differences  in  the  average  score  profiles  of 
the  four  groups. 

>  Establishing  procedures  for  classifying  boating  trajectories 
on  the  basis  of  their  score  on  a  set  of  attributes. 


> Standardized  Canonical  Discriminant  Functions: 

Zy  —  0.047xj  +  1.009x2  +  0.585x3  +0.093x4 
Z2  =  0.454x1  -0.584x2  +0.469x3  +0.671x4 
Z3  =  -0.744XJ  +0.202x2  +0.732x3  +0.174x4 

>Fisher’s  Linear  Classification  Functions: 

Canoe  = -6.941-5.223£-05X,  +  3.301X,  +  0.125X,  +  5.894£-04X4 
Kayak  =  -5.629  - 1 .5 16 E-  05X,  +  3.25  \X2  +  0.1 03X,  +3.209£  -  04X„ 
Motorboat  =  -27. 1 6  -  6.223 E  -  05X,  +  9.891X2  +  0. 19 1X3  +  7.6 17£  -  05X4 
Sailboat  = -13.74  +  3.057.E-05X,  +  5.156X2  +  0.179X,  +1.498£-03X4 


xl=Total  distance,  x2=Mean  speed,  x3=mean  turning  angle,  x4=farthest 
distance  to  shore. 


Predicted  Group 

Membership 

Total 

TYPE 

canoe 

kayak 

motorboat 

sailboat 

Original 

Count 

canoe 

15 

2 

0 

0 

17 

kayak 

3 

18 

0 

0 

21 

motorboat 

1 

0 

8 

1 

10 

sailboat 

0 

0 

0 

47 

47 

% 

canoe 

88.2 

11.8 

.0 

.0 

100.0 

kayak 

14.3 

85.7 

.0 

.0 

100.0 

motorboat 

10.0 

.0 

80.0 

10.0 

100.0 

sailboat 

.0 

.0 

.0 

100.0 

100.0 

92.6%  of  original  grouped  cases  correctly  classified. 


Key  results 


>  Stepwise  discriminant  analyses  for  the  coastal  and  river 
study  areas  retain  the  same  independent  variables 

>Mean  Speed,  Mean  Turning  Angle,  Farthest  distance  to 
shore 

> Geography  does  not  influence  boating  trajectories  much 

>Mean  Speed,  Mean  Turning  Angle,  Farthest  distance  to 
shore  are  significantly  different  among  different  groups 


P(Kayak)=0.733 

P(Canoe)=0.240 


P(Motorboat)=0.899 
P(Sailboat)  =0.101 
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